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Clusters of galaxies are embedded in halos of optically thin, gravitationally stratified, 
weakly magnetized plasma at the system’s virial temperature. Due to radiative cooling 
and anisotropic heat conduction, such intracluster medium (ICM) is subject to local in¬ 
stabilities, which are combinations of the thermal, magnetothermal and heat-flux-driven 
buoyancy instabilities. If the ICM rotates significantly, its stability properties are sub¬ 
stantially modified and, in particular, also the magnetorotational instability (MRI) can 
play an important role. We study simple models of rotating cool-core clusters and we 
demonstrate that the MRI can be the dominant instability over significant portions of 
the clusters, with possible implications for the dynamics and evolution of the cool cores. 
Our results give further motivation for measuring the rotation of the ICM with future 
X-ray missions such as ASTRO-H and ATHENA. 

PACS codes: 


1. Introduction 

Galaxy clusters are embedded in gaseous halos at the system’s virial temperature 
(10^ — 10® K). Such intracluster medium (ICM) is a dilute, weakly magnetized plasma, 
stratified in the cluster gravitational potential. The question of the local stability of the 
ICM is a fundamental piece in the puzzle of the evolution of clusters and in particu¬ 
lar of cool-core clusters, in which the cooling time of the plasma in the central regions 
is shorter than the cluster age. The results of linear-stability analysis and magneto¬ 
hydrodynamics (MHD) simulations indicate that, due to the joint effect of radiative 
cooling and anisotropic heat conduction, the ICM is subject to several local instabil¬ 
ities, which are combinatio ns of the the rmal instability (TI; iFie 13 Il965l) . the magne- 
tother mal instability (MTI; I BalbusI [200011 and the heat-flux-driven buoyancy instability 
(HBI; lOnataert 2008h . Cosmological hydrodynamic simulations suggest that the X-ray 
emitting halos of clusters, though mainly pressure supported, might rotate significantly 
( Fang fit al. 20091 : Fan eJ, alll2012h . with a contributio n of rotation especially important 
in the cool cores of relaxed clusters ( Nagai et al. I l201,3l . Observationally, the most direct 
signature of ICM rotation would be the spatially resolved measure of shifted X-ray emis¬ 
sion line centroids. Unfortunately, given the relatively poor spectral resolution (> 100 eV) 
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of currently available X-ray instruments, so far there are no meaningful constraints on ro¬ 
tation based on such a measure, which requ ires the spectral capab i lities of the up-coming 
X-ray calorimeters on board of ASTRO-H ( Takahashi et a/.ll2014l : iKitavama et al\\201^ 
and ATHENA ( Ettori et al. 2013h . Spatially unresolved rotation contributes, together 
with the turbulence, to the emissi on line broadening: the present observation al measures 
of the X-ray emission line widths ( Sanders fc Fabian 20131 : Pinto et all 2015h . even when 
combined with m easures of the flatteni ng of the X-ray iso photes, leave ample room for ro¬ 
tational motions ( Bianconi et a/.ll2013L hereaft er BEN13 ). If the I CM rotates, its stability 
properties are substantially modified ( Ba.lbus 2001 : Nipotil 2010ll and, in particular ^ also 
the magnetorotational ins tability fMRI: iBalbus fc Hawlm; 1^91 ; see als o VelikhovI 19591 
and Chandrasekhar! 1960f l can play an important role ( Niooti fc Postil 2014 , hereafter 
NPm 

The linear evolution of local axisymmetric perturbations in a stratified, rotating, ra- 
diatively c ooling, weakly magne tized plas ma is d e termin ed by the dispersion relation 
derived in iNipoti fc Posti ( 2(1131 hereafter NP13h . NP14 studied the nature of the in¬ 
stabilities of such a plasma in two illustrative cases, generically representative of the 
physical conditions in galactic and cluster atmospheres. Here we focus on the cool cores 
of galaxy clusters and we disc uss the s tability properties of rotating cool-core models 
built on the basis of the work of lBENl^ The present study extends the results of lNPl4 
as here we study the nature and linear growth rate of the plasma instabilities through¬ 
out observationally motivated cool-core cluster m odels. In this work we limit ourselves 
to axisymmetric perturbations, but the results of jNPlJ suggest that the study of non- 
axisymmetric disturbances should lead to similar conclusions. Though we consider global 
cluster models, we study only the stability against local perturbations, i.e. disturbances 
with sizes much smaller than the characteristic scale-lengths of the system. In addition 
to the local instabili ties found in the pres ent study, the ICM could be subject to global 
unstable modes (e.g. Latter fc Kund 2012 ). which by construction would elude our anal¬ 
ysis. 

The paper is organized as follows. The governing MHD equations are given in Sectionj^l 
the unperturbed models are presented in Section [3] and the dispersion relation obtained 
from the linear-stability analysis is given in Section |4l Section [5] describes the properties 
of the local instabilities occurring in the cluster models and Section [ 6 ] concludes. 


2. Governing equations 

A stratified, rotating, magnetized, dilute plasma in the presence of thermal conduction 
and radiative cooling is governed by the following MHD equations: 
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with the additional condition V • B = 0. Here p, p, T, v and B are, respectively, the 
density, pressure, temperature, velocity and magnetic field of the fluid, $ is the external 
gravitational potential (we neglect self-gravity), 7 = 5/3 is the adiabatic index, H is the 
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heating rate per unit volume, C = CiT, p) is the radiative energy loss per unit mass of 
fluid, and 


Q = - 


xB (B-v)r 

52 


(2.5) 


is the conductive heat flux, where y = kT ®/2 jg ISpitzeil (jl962ll electron conductivity 


with K, ~ 1.84 X 10“^(lnA) ^ ergs“^ cm“^ and In A is the Coulomb logarithm. 

Neglecting the weak temperature and density dependence of In A, here we fix In A = 30, 
so K is a constant and x — x{T) oc As we consider rotating fluids, we work 

in cylindrical coordinates {R, (j), z), where i? = 0 is the rotation axis. The governing 
equations (12.1112.41) are given explicitly in cylindrical coordinates in Appendix [A) 

The expression (12.51) for Q accounts for the fact that in a dilute magnetized plasma 
heat is significantly transported by electrons only along the magnetic field lines. For 
simplicity, even if the medium is magnetized, we have assumed that the pressure is 
isotropic, neglecting the fact that momentum transport is anisotropic in the presence 
of a magnetic field (i.e. the so-called Braginskii viscosity; Brasfinsl^ll965l ). Though this 
approximation is not rigorously justified, we adopt it in the working hypothesis that 
anisotropic pressure is not the crucial factor in determining the dominant instabilities 
of a rotating ICM. The importance of the effects of the Braginskii viscos ity on plasma 
stability is debated ( Kund 2011 : Latter fc KunzlEoi 2 l : IParrish et al. 2012 ). so the limits 
of this approximation must be taken into account when interpreting the results of the 
present work. 

The term H in equation (12.41) represents an unspecified source of heating that is ex¬ 
pected to counteract radiative cooling and prevent a cooling catastrophe. There is general 
cons ensus that heat con duction cannot be entirely responsible for halting cooling flows 
(e.g. Parrish et al ][ 2 ^ and that other sources of heating, such as feedback from the 
central active galactic nucleus (AGN), must be at work to balance cooling in a time- 
averaged sense. In the absence of physically motivated and analytically tractable models 
of such a heating term, here we just consider a toy model in which, f or a given clus¬ 
ter m odel, 7^ is a fixed position-dependent function (for a discussion see iMcCourt et al. 


2OI2I) : the main results of this work are independent of this approximation, though very 
crude. 


3. Cool-core cluster models 

We consider here three rotating cool-core cluster models (n amed Cl, C2 and C3), which 
are slightly modified versions of those presented in BEN13I (indicated as CCl, CC 2 and 
CC3 in that paper). In all these models, which represent massive galaxy clusters, the total 
gravitational potenti al $ is given by a sphericaQ Navarro-Frenk-White (NFW) model 
( Navarro et a/.l[l9^ with virial mass M 200 = scale radius r^ = 519kpc and 

virial radius r 2 oo = 2066 kpc. The plasma rotates differentially with azimuthal velocity 
ro 0 = VLR, where = VL{R) is the a ngular v elocity and the rotation law is different in 
the three models (see section 2.3 of BEN13I ). The other velocity components are null. 
Assuming that the ICM magnetic field is dynamically unimportant, we construct the 
unperturbed pressurepo = PoiR, z), density po = po{R, z) and temperature Tq = Tq{R, z) 
fields as axisymmetric stationary solutions of equations ( 12 . 1112 . 21 ) with B = 0. As 11 does 
not depend on z, the gas distribution is barotropic {po, po and Tq are stratified on 


f In a more realistic model the total gravitational potential would be axisymmetric or triaxial. 
Though idealized, a model with spherical potential is an interesting limiting case because the 
flattening of the ICM distribution is entirely due to rotation. 
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log (R/kpc) 

Figure 1. ICM electron number density profiles in the equatorial plane of the rotating 
cool-core cluster models Cl, C2 and C3 (see Section [3]|. 



log (R/kpc) 


Figure 2. ICM temperature profiles in the equatorial plane of the same models as in Fig. [T] 
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log (R/kpc) 

Figure 3. ICM rotation speed as a function of radius for the same models as in Fig. [T] 
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Figure 4. ICM angular velocity as a function of radius for the same models as in Fig. [T] 
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log (R/kpc) 

Figure 5. ICM magnetic field modulus as a function of radius in the equatorial plane of the 

same models as in Fig. [T] 


the surfaces of constant effective potential). In particular, the plasma is modeled as an 
equilibrium two-component composite polytrope with outer polytropic index %ut = 1.14 
and inner polytropic index yin < 1, so the ICM temperature decreases outwards in the 
outer regions and increases outwards in the inner regions (cool core). The parameters of 
m odels Cl , C2 and C3 are, respectively, the same as those of models CCl, CC2 and CCS 
of BENISL with the only exception of the inner polytropic index yin, which is yin = 0.49 
in model C l, yin = 0.43 in model C2, yin = 0.56 in model C3 (these values, smaller than 


BENlSl . lead to higher central densities and lower central temperatures). The plasma 


temperature is in the range 2 — 7k eV and the size of the cool core is about 140 kpc in 
the equatorial plane. As in IBENISI we assume that the ICM metallicity is everywhere 
3/10 of the solar metallicity. Figs. [TJUlplot the radial profiles, in the equatorial plane, 
of electron number density Ueo, temperature Tq, rotation speed vo^j, and angular velocity 
n. As discussed in IbENISI the above models of rotating clusters are realistic, in the 
sense that they are consistent with the current measurements of the ellipticity of the 
X-ray isophotes and of the width of the X-ray emission lines. Additional constraints 
on the rotation of the ICM derive from the so-called hydrostatic mass bias, that is the 
discrepancy between cluster mass estimates from gravitational lensing and those based 
on the assumption that the ICM is in hydrostatic equilibrium. In Appendix [B] we show 
that models Cl, C2 and C3 are realistic also in this respect. 

Given a cluster model, we assume that the heating term T-L = 'H{R, z) appearing 
in the energy equation ()2.4I) is such that H — poC{To,po) — V • Qo =0, where Qo 
is the unperturbed conductive heat flux. We remark that in our models Ti is not a 
function of the hydrodynamic variables, but depends explicitly on the coordinates R and 
z. To complete our cluster models we must specify the properties of the unperturbed 
magnetic field Bq, which, though not relevant to the equilibrium configuration (because 
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P = Sttpq/Bq » 1), influences the stability of the ICM. We do not attempt a full 
characterization of the magnetic field throughout the global cluster model. For simplicity, 
In this work we focus only on the equatorial plane z = 0, and we specify locally Bg as a 
function of cyl indrical radius R. The stationary unperturbed magnetic field must satisfy 
Ferraro! (Il937r) isorotation law Bg • = 0, so Bqr = 0 in the considered barotropic 


distributions. For each cluster model, we consider two cases: one in which the only non- 
vanlshlng component of the unperturbed magnetic field Is Bq^ (so Bq^ = 0) and another 
in which both Bg^ and Bq^ are non-vanishing {Bq^/Bq^ = 10). In all cases, given the 
unperturbed pressure distribution pg, the unperturbed magnetic field modulus Is fixed 
by imposing a constant valu e of P such tha t at the centre Bq = 10/iG (consistent with 


observational estimates; e.g. Brhggen 201.3l l. l.e. P = 29.3 for model Cl, P = 26.8 for 


model C2 and P = 89 for model C3. The profiles of the magnetic field modulus for the 
three models are shown In Fig. [S] We note that in the equatorial plane the unperturbed 
conductive heat flux Qg vanishes, because at z = 0 the magnetic field lines are isothermal 
(the only non vanishing component of the temperature gradient is the radial one, which, 
combined with Bg^j = 0, implies Qg = 0), so the condition of energy balance is simply 
^ = Po^(To, po). 


4. Linear stability analysis 

The stability analysis is performed by linearizing the system (I2.1II2.4D with local ax- 
isymmetric Eulerian perturbations of the form with |/| |/g|, where 

fo is the unperturbed quantity, cu is the perturbation frequency, and kjt and kjz are, re¬ 
spectively, the radial and vertical components of the perturbation wave-vector k. In the 
short wave- length and low frequency approximation, the resulting dispersion relation for 
n = —iw is (NP13) 


+ oJdTT-'^ + [^Bv T ^rot + + [(Wfot + 2cij^)a;d -I- w^Wc.mag] 


+^A ( + ^BV + ^rot “ ^ + 


W, 


c,<p 


n 




wl- 


' ^rot 






I ^ 

“T CdA^c,mag 


= 0 , 


(4.1) 


where the quantities indicated with subscripted w are characteristic frequencies related to 
local properties of the unperturbed distribution and of the wave-vector. Specifically, wg = 
Wth + Wc,a, where Wth is the TI frequency and uJc,a is the anisotropic thermal conduction 
frequency, wbv is the Brunt-Vaisala frequency, Wrot is the angular momentum gradient 
frequency, uj\ is the Alfven frequency, and Wc^mag and are two other frequencies 
associated with th ermal conduction mediated by the magnetic field (all the definitions 


are give n in NP13h . The radiative energy loss per unit mass of fluid £ (pg, Tg), appearing in 


cut vi (see NP13|), is computed u sing the collisional ionization equilibrium cooling function 


of Sutherland & Dooita (1993). In the disoersion relation (ED there are no terms related 
to the heating source 77 appearing in equation (12.41) , because we are considering Eulerian 
perturbations at a fixed point in space and 77 is a function of position, but not of the 
hydrodynamic variables. For fixed wave-vector k, given a cool-core cluster model (see 
Section ID, the coefficients of equation ED are fully determined, so the problem of 
the local linear stability of the plasma is reduced to finding the roots of a polynomial 
equation. 
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5. Results 

For each of the considered cluster models, at different radii R in the equatorial plane, for 
given wave-vector k we numerically computed the solutions of the disp ersion relation 63) 
using the IDL (Interactive Data Language) routine fz_roots (as in INP141 ). Depending 
on the signs of their real (R„) and imaginary (I„) parts, the solutions are classified as 
stable (R„ ^ 0), overstable (R„ > 0, I„ 7 ^ 0) or monotonically unstable (R„ > 0, I„ = 0) 
modes. As overstable linear disturbance s are unlikely to enter the non-linear regime 
( Malagoli et aI.IUqS? : Binnev et al. 2009| l. here we focus only on the more interesting 
monotonically unstable linear modes. For each model we compute, as a function of radius 
R, the growth rate ninst of the fastest-growing monotonically unstable mode. At each 
radius, riinst is calculated considering modes with wave-number k = + in the 

range IQ < kh < 1000, where h = (|Vp|/p)“^ is the pressure scale-height, which depends 
on R\ this guarantees that the perturbation wave-length is short with respect to the 
characteristic scales of the system at all radii. The radial-to-vertical wave-number ratio 
kn/kz is allowed t o vary in t he ran ge —25 ^ ka/kz ^ 25. 

As discussed in lNP13l and InP 14 in general the instabilities of rotating, weakly mag¬ 
netized, radiatively cooling plasmas arise from combinations of TI, MTI, HBI and MRI. 
In order to understand the physical consequences of t he inst ability it is important to 
determine the nature of the dominant unstable modes. iNPldl have shown that, though 
the different modes are entangled, it is possible to identify the different branches of the 
solutions in terms of combinations of well-known modes (TI, MRI, MTI, HBI, rotation, 
buoyancy and Alfven modes). In practice, this is done by inspecting the behaviour and 
the growth rates of the individual branches of the solutions as functions of the wave- 
number and comparing the results with those obtained for simpler configurations (for 
instance in the absence of radiative cooling, rotation or magnetic field). Applying this 
technique to the cluster models here considered, we have been able to identify the nature 
of the dominant instabilities at different radii. 

It is useful to compare the maximum instability growth rate riinst with the characteristic 
rates of radiative cooling and of the MRI. Therefore, we define the cooling rate 


_ 7 - 1 /iTOp£(To,po) 

'^cool — 


(5.1) 


7 ^bTo 

which for each unperturbed model is a function of position (here /i = 0.59 is the mean 
mass per particle in units of the proton mass rup). The MRI growth rate timri is the max¬ 
imum gr owth rate of monotonica lly unstable modes that are solutions of the dispersion 
relation ( Balbus fc Hawlev 199lh 


+ [‘^BV ' 


' ^rot 


■ 20;! 


■uj\ 




^BV ' 


A^ot-4f22^)=0, 


(5.2) 


which is obtained from the more general dispersion relation 63 in the limit of absence 
of radiative cooling and thermal conduction (wd = 0, Wc,mag = 0, = 0). As well 

known, the TI growth rate |wth| is of the order of, but typically somewhat higher than 
the cooli ng rate tirn r,], which is at the basis of the existence of the TI in astrophysical 
plasmas (IField 1965 ). The MRI gro wth rate timri is proportional to the plasma angular 
velocity ( Balbus fc Hawle\^ 1 1991 ). 

As mentioned in Section |31 for each cluster model, we consider two cases: one in which 
Bo 4> = 0 and another in which Bq^/Bqz = 10 at all radii. The only difference between 
these two families of model is the value of the ratio Bq^/Bqz- all the other properties are 
the same, including the radial profile of the magnetic field modulus, which is shown in 

Fig. El 
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Figure 6. Maximum instability growth rate riinst, cooling rate ricooi, and MRI growth rate timri 
as functions of radius in the equatorial plane of the same models as in Fig.[T] Here = Bqr — 0 
and Boz ^ 0. At R > 80kpc there are no monotonically unstable modes, at least in the explored 
wave-vector space. 



log (R/kpc) 


Figure 7. Instability (tinst), cooling (tcooi) and MRI (tMRi) timescales normalized to the 
free-fall time tg as functions of radius in the equatorial plane of the same models as in Fig. 
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5.1. Models with vanishing azimuthal magnetic field (Bq^ = Oj 

Fig. [6] plots, as a function of radius, the maximum growth rate n-mst in the equatorial 
plane of models Cl, C2 and C3 when Bq^ = 0. For comparison, the cooling (ucooi) 
and MRI (nMRi) rates are also plotted for the same models. All models are unstable 
in their cool core, while no monotonically unstable modes are found at i? > 80 kpc. 
Model Cl is unstable for MRI in the inner ^ 10kpc (instability timescale ~ 10®yr), 
while at larger radii the system is unstable for TI (instability timescale ~ 10® yr, of the 
order of the cooling time). Models C2 and C3 are not unstable for MRI, but present 
TI modes throughout the cool core with instability timescales between 5 x 10® and 
5 X10® yr. Together with the absolute timescales, it is also useful to consider the timescales 
normalized to the local dynamical (“free-fall”) time te = [2r/g(r)]^/®, where g = dd)/dr 
is the modulus of the gravitational force and r is the radial spherical coordinate. Fig. [7] 
plots, as functions of radius, the instability (tmst = ^inL)) cooling (tcooi = 

MRI (^mri = ^MRi) timescales normalized to tg. It is apparent that, when present, the 
MRI grows on timescales as short as 5tg, while the TI timescale is typically 10 — lOOtg. 

Though all models satisfy at all radii the necessary condition to have MRI (Rqz 0 
and dn/di? < 0; see Fig. 0]), MRI modes are present only in the inner core of model 
Cl. We recall that the MRI occurs only for wave-numbers smaller than a critical wave- 
number /cmri, which, with the exception of the centre of model Cl, is outside the explored 
wave-number interval (in other words, short wave-length perturbations are not unstable 
for MRI in these cases). Neither the HBI nor the MTI contributes to the unstable modes: 
these instabilities are related to anisotropic thermal conduction and can occur when the 
temperature either increases (HBI) or decreases (MTI) outwards. In the cool core the 
temperature increases outwards, but the magnetic field lines are isothermal so no HBI is 
expected. Out of the cool core the temperature decreases outwards, so in principle the 
MTI might occur, but in fact no MTI modes are found in the explored wave-vector space. 
Similarly to the MRI, the MT I occurs only for wave-numbers smaller than a critical wave- 
number A:mti < Balbus 2000ll . which is out of the explored wave-number interval (as it 
happens for the MRI, short wave-length disturbances are not unstable for MTI in these 
cases). 


5.2. Models with dominant azimuthal magnetic field (Bqi^/Bq^ = 10^ 

Figs. [HE] plot for models Cl, C2 and C3 the same quantities as Figs.|6]|7l but under the 
assumption that the azimuthal magnetic field component is dominant {B 04 ,/Bqz = 10). In 
this case all models present monotonically unstable modes out to large radii (~ Mpc). 
At all radii the unstable modes of model Cl have growth rate much higher than the 
cooling rate: the MRI dominates throughout the cool core, while in the cluster outskirts 
the unstable modes are driven by the MTI. Model C2 presents TI modes in the inner 
parts of the cool core, MRI modes in the intermediate radial range (10 — 100 kpc) and 
MTI modes in the outskirts. Model C3, which has the lowest rotation speed at all radii 
(Fig. |3|) is stable against MRI even when Bq^ is dominant: its unstable modes are driven 
by TI in the cool core and by MTI in the outer regions where the temperature gradient 
is negative. In the explored models, both the MRI and the MTI modes have remarkably 
short timescales, of the order of a few local dynamical times. In absolute terms the MRI 
modes are fastest with growth times as short as ^ 10® yr in the core. 

Comparing Figs.|HE|with Figs.|n]|7|it is apparent that the instabilities are stronger and 
more widespread when the magnetic field is dominated by the azimuthal component. The 
driving factor in determining these differences is not the presence of an azimuthal mag¬ 
netic field component, but the fact that when Bq^ is dominant, the vertical component of 
the magnetic field is weaker (at fixed Bq). As pointed out above, the MRI and the MTI 
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Figure 8. Maximum instability growth rate Uinst, cooling rate ricooi, and MRI growth rate timri 
as functions of radius in the eqnatorial plane of the same models as in Fig. [T] Here Bor = 0, 
Boz 7^ 0 and Bq^/Boz = 10. 



log (R/kpc) 

Figure 9. Instability (finst), cooling (fcooi) and MRI (tiuRi) timescales normalized to the 
free-fall time ts as fnnctions of radius in the equatorial plane of the same models as in Fig. |8l 
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occur only for wave-numbers smaller than critical values (/cmri and kuTi, respectively). 
Both /cmri and ^mti are inversely proportional to the magn etic field component coupled 
with the wave-vector ( Balbus fc Hawlev 1991 ; BalbusI 2000ll , which in the present case is 
just i?oz- When the azimuthal magnetic field is dominant, the vertical magnetic field is 
sufficiently weak to trigger the MRI throughout the cores of models Cl and C2, and the 
MTI in the outer regions of all models. It is also the case that ^mri increases for increas¬ 
ing angular velocity gradient, which explains why the MRI manifests itself in models 
Cl and C2, which have stronger angular velocity gradients than model C3 (see Fig. S]). 
When present, both the MRI and the MTI have growth rates much higher than the local 
cooling rate: the fastest instability is the MRI, which acts on timescales 10® — 10® yr. 


6. Conclusions 

Building upon the results of NPI3L NP14 and BEN13 . in this paper we have studied 
local instabilities arising from axisymmetric perturbations in rotating cool-core clusters. 
For the sake of simplicity we limited ourselves to somewhat idealized cluster models. 
Notwithstanding these limitations, our calculations should catch the essential properties 
of more realistic cluster models that could be built in the near future, especially when 
more observational constraints on the ICM rotation are available, thanks to the next- 
generation X-ray missions such as ASTRO-H and ATHENA. 

The results of the linear stability analysis of the rotating cool-core cluster models here 
presented show that, depending on the ICM properties, TI, MTI and MRI can occur 
at different radii in the equatorial plane (where, by construction, our models exclude 
the HBI, which however can play a role in more general configurations). What are the 
astrophysical consequences of these instabilities? 

In principle, monotonically unstable TI modes can lead to local condensation of cold 
gas out of the hot plasma and therefore to a multiphase ICM. However, given that the 
growth rate of TIs is of the order of the cooling rate, the evolution of thermally unstable 
modes depend on the details of the balance of heating and cooling in the unperturbed 
fluid. While heat conduction always tend to damp thermal perturbations, other heating 
mechanisms such as AGN feedb ack have more complex effects and can in som e cases even 
favour cooling instabilities le.g. ICiotti & Ostrikeiil2007t ICasoari et aLll2012ll . Therefore, 
it is not possible to draw robust conclusions on condensational modes based on the models 
here considered, in which we have adopted a very simple ad hoc heating source to balance 
radiat ive losses in the background plasma. Multiphase plasma is detected in cluster cool 
cores (Heckman I98lll . but the origin of the colder medium is still debated. The observed 
correlation betwe en presence of multiph ase gas and short cooling time (low entropy) in 
cluster cores (e.gJVoi£.^,Donahue j201® does not necessarily imply that the cold gas is 
produced by TI ( Niooti fc BinnevI 20041 1. 

As far as the MTI is concerned, our results confirm previous findings that the MTI 
can be the dominant instability out of the cluster core. As in non-rotating ICM models, 
the MTI can affect the outer regions of the cluster (where the temperature decreases 
outwards) on times scales as short as a few local dynamica l time s, potentially driving 
strong turbulence in the non-linear regime ( McCourt et al.lEoil l. We recall that our 
models do not account self-consistently for anisotropic momentum transport, which in 
principle can affect the behaviour of the MTI. However, it has been shown that, at least 
in non-rotat ing systems, the effec ts of the Braginskii viscosity, which can be significan t 
for the HBI ( Latter fc Kunz 2012 1. are almost negligible for the MTI ( Kunz et al. 2012 1. 

The most interesting instabilities among those found in this work are the MRI modes 
in the cool cores. Figs. [BE] show that the MRI rates are remarkably high, with associated 
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timescales 10® — 10® yr (a few local dynamical times): this suggest that the MRI can have 
an important impact on cool cores. The MRI timescales are substantially shorter than the 
cooling time, so this conclusion is robust against uncertainties on the balance of heating 
and cooling in the core. The Braginskii viscosity, which is neglected in o ur calculations, 
could even enhanc e the growth rates of these dominant MRI mo des (lOuataert et all 
2 OO 2 I : [BalbusI 12004 1. Studies of weakly magnetized accretion discs ( Balbus et all 199 '^ 


have shown that the MRI is expected to drive turbulence, with associated energy and 
angular momentum transport. These MRI-driven dissipative processes can contribute 
substantially to the dynamics and evolution of cluster cool cores, if they rotate signif¬ 
icantly. The MRI could be the engine of turbulent heating, which is believed to be an 
efficient mechanism to halt cooling flows in galaxy clusters ( Zhuravleva et al. 20I4h . A 
natural development of the present work will be a quantitative estimate of the effect of 
the MRI on the energy balance and accretion rate of the plasma in cluster cool cores. 
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thank the two anonymous referees for useful comments. CN has been supported by PRIN 
MIUR 2010-2011, project “The Chemical and Dynamical Evolution of the Milky Way 
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Appendix A. Governing eqnations in cylindrical coordinates 

In cylindrical coordinates the governing equations (12.1112.41) read 


dp 1 dRpvn 1 dpvtj) dpvz 

dt~^R dR R d4> dz 
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where we have assumed that the gravitational potential is axisymmetric and we have 
used V • B = 0 in writing the three components (equations lA 5IIA 71) of the induction 
equation dSSl). 
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Figure 10. Hydrostatic (Mtherm) to true (Mtme) mass ratio within a sphere of radius r for the 

same models as in Fig. [T] 


Appendix B. Hydrostatic mass bias 

The estimate of the mass of galaxy clusters is mainly based on X-ray and gravitational 
leasing data. Mass estimates from X-rays are obtained in general under the assumption 
that the gas is in hydrostatic equilibrium, while mass estimates from gravitational leasing 
do not depend on the kinematics of gas. Deviations from hydrostatic equilibrium (due to 
rotation, bulk motions or turbulence) produce a discrepancy between the mass estimated 
from X-ray data and the mass estimated from gravitational leasing. Such a discrepancy, 
which is usually called hydrostatic mass bias, can be quantified by measuring within a 
sphere of radius r the ratio Mtherm (u)/Mtrue(u), where Mtrue is the true mass and Mtherm 
is the m ass computed as suming that gravity is balanced only by the thermal-pressure 
gradient ( Lau et al 

The ratio Mtherm/Mtrue is estimated in hydrodynamic cosmological simulations to 
be around 0.8 — 0. 9, with a slight variation with radius out to r 2 oo and a scatter of 


When X-ra y simulated 


about 20 per cent (iMeneghetti et al. 2010[ iNelson et al 
observations are used to mimic typical X-ray exposures, iRasia et ol\ (|2012l l show that 
temperature inhomogeneities could affect further the hydrostatic mass reconstruction 
by reducing Mtherm/Mtrue to 0.65 — 0.75. On the observational side, we cannot probe 
the “true” mass but only a reasonable proxy for it as pr ovided by gravitational lens - 
ing esti mates. Programs like Weighing the Giants (WtG, Ivon der Linden et al.ll20i3 l. 
CLASH ( Dona hue et al JHoii, and the Canadian Cluster Comparison Project (CCCP, 
Hoekstra et aLll2015[l obtain values of the hydrostatic-to-lensing mass ratio in the range 


of 0.7 — 0.8, although differences up to 40 per cent in either the weak leasing or X-ray 
mass measurement among different research groups are still measured (ISereno &: Ettori 
201511 . 

For models of rotating cluster such as Cl, C2 and C3 (see Section [3]) Mtherm < Mtme, 
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because of the presence of rotational support. Models in which rotational support is more 
important are characterized by lower values of Mtherm/Mtme' for a rotating model to be 
realistic, it must have Mtherm/Mtme higher than the aforementioned observed values. It is 
therefore interesting to estimate the hydrostatic mass bias for models Cl, C2 and C3. In 
all three cases Mtrue(c) is just given by the spherically symmetric N FW mass distrib ution 
responsible for the total cluster potential (see Section [3]). Following iLau et d. ( 2013h . the 
hydrostatic mass profile can be computed as 


.Mtherm(r') 


d (po) 

G (po) dr 


(Bl) 


where the symbol (■ • •) indicates an average over angles at fixed spherical radial coordi¬ 
nate r. We construct (po) and (po) by interpolating pq{R,Z) and po{R, z) in spherical 
coordinates (r, 9) and then computing the average over 0, at fixed r. The profiles of 
Mtherm/Mtrue for models Cl, C2 and C3 are shown in Fig. [TUI The minimum values of 
A^therm/foftrue are 0.76, 0.66 and 0.94, for models Cl, C2 and C3, respectively. However, 
these minima occur within the cool-core region (r < 100 kpc), in which the hydrostatic 
mass bias is not easily measured observationally. The minimum values of Mtherm/foftrue 
in the radial range Tg ^ r ^ r 2 oo (he. 0.5 < r/Mpc < 2), which is representative of the 
region typically probed by weak gravitational leasing, are 0.88 (Cl), 0.90 (C2) and 0.98 
(C3), consistent with the observational limits. 
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